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Abstract 

Using the Post-Minkowskian formalism and considering rotation as a 
perturbation, we compute an approximate interior solution for a sta- 
O ■ tionary perfect fluid with constant density and axial symmetry. A 

suitable change of coordinates allows this metric to be matched to the 
exterior metric to a particle with a pole-dipole-quadrupole structure, 
relating the parameters of both. 

o : 
<n : 

7— I . 

1. Introduction. 

. 

In the study of multipole particles in General Relativity, certain approximate solutions 
■ to the gravitational field outside such objects have been proposed. The information about 
their structure is collected in a series of parameters, which are interpreted as the multipole 
moments of the particle. Here, we use as a starting point the metric outside a particle 
with pole-dipole-quadrupole structure (PDQ particle), with stationary and axisymmetric 
moments [1]. Our aim is to find an approximate interior metric that will match it on a 
given surface. It is demanded that they should match in the sense of Lichnerowicz [2] : the 
metric and its first derivatives must be continuous on the matching surface. For simplicity's 
sake we shall assume that the inside is made up of a homogeneous perfect fluid in rigid 
and slow rotation (this term will be precised below). The interior solution will be obtained 
by solving the Einstein equations using the Post-Minkowskian perturbation method [3]. 
From study of the motion equations of the material medium it is possible to extract precise 
information for obtaining the equation of the matching surface [4]. Finally, the matching 
conditions of both metrics will be translated into relationships between the parameters of 
the interior metric (density and angular speed of rotation), the mean radius of the surface, 
and the parameters of the exterior metric (multipole moments). 

2. The exterior metric. 

In a system of harmonic coordinates X a such that X° = t and <p are adapted to the 
Killings vectors (R, O, and <p are the spherical coordinates associated to the harmonic 
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ones), the exterior metric to an isolated PDQ particle with a stationary and axisymmetric 
structure is written as follows (see ref. 1): 
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where 



E l = 5\ , N l = — , M l = e i jk N j E k 



(2-2) 



Pz(cosG) is the Legendre polynomial of order / and the parameters m, S and B, are 
respectively the mass, the angular momentum and the mass quadrupole moment of the 
particle. 

Note that the terms containing the different moments have a specific dependence on 
the powers of R~ x . More precisely, it may be seen that only the following factors appear: 
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and their powers and products. 

Finally, it is worth mentioning that (2.1) coincides with the expansions, in certain 
harmonic coordinates, of Schwarzschild, Weyl and Kerr metrics, if the parameters m, S 
and B take the values: 



Gm = a GS = GB = (Schwarzschild) 
Gm = -ao GS = GB = |ag - a 2 (Weyl) 
Gm = a GS = —a a GB = —a a 2 (Kerr) 
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(2-4) 



3. Newtonian theory. 



The Post-Minkowskian approximation allows one to obtain the interior metric in the 
form of an expansion in powers of G, which we shall match, at each iteration, with the 
terms of the same order of the exterior solution. Since the results obtained at first-order 
must necessarily coincide with those afforded by Newtonian mechanics, it makes sense to 
first study the problem posed within this framework. 

It should be recalled that study of a self-gravitating fluid is already a non-linear 
problem in Newtonian mechanics. On one hand, the field equation: 

Aip G = A-KGp (3.1) 

allows one to obtain the gravitational potential, ipa, once we know the mass distribution 
p. On the other hand, Euler motion equations: 

d t p + d k (pv k )=0 

(3.2) 

d t {pv l ) + d k (pv k v* + a kl ) = -pd^ G 



(where v l is velocity and a %k is the stress tensor of the fluid) afford information about the 
evolution of the source of the gravitational field, if we have prior knowledge of the potential 
ipG inside it. 

Let us assume that the fluid is perfect (a lk = p5 lk ) : that density is constant in time 
{dtp = 0) and that the fluid rotates in a rigid manner with a constant angular velocity lo 1 
{v l = e l X k ). Under these conditions, equation (3.2) can be reduced to (we choose the 
X 3 axis in the direction of 

d^p = , dip = -pd^ (3.3) 

with 

tj, = i(;g - -wV sin 2 6 (3.4) 

The second equation in (3.3) is a constraint on p, p and ip. If one considers this as an 
equation for p, from its integrability condition 



d [l pd j] il) = Q (3.5) 

it may be deduced that p must be a function of ip, p = x(V0 ( or the other way round); 
and its solution is: 

p = p - I X(A) d\ (3.6) 

Jipo 

The surface E of the fluid is determined by the equilibrium condition p |s= 0, which leads 

p = [ S X (A) d\ (3.7) 
J4> 
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where is the constant value of tp on E (hence, ip(R, 0) = V>e is the implicit equation 
of S). Moreover, since p and p are functions of ip, they are constant on the surfaces 
ip(R, 0) = constant. 

If an equation of state p = f(p) is given as an input, (3.3) allows one to obtain a 
relation p = x(V0> which; substituted in (3.1), leads to a generally non-linear equation for 
the potential: 

A^ G =4^x0/0 (3-8) 
For example, for a polytropic-like state equation 

p = a p 1+ ~ (a = constant , 7 > 0) (3.9) 



one obtains 
and 

A^ G = 4tvG C - (3.11) 

which is not linear if 7 7^ 1 [5] . 

The case of p = constant, which is not included in the previous example, is not even 
free of non-linearities. We have already mentioned that ip(R, 0) = i/j'E defines the surface 
E. However, to determine it we must first know the potential tjjQ, which in turns depends 
on the shape of the surface 

V>g = -G I - 9 - d 3 X' (3.12) 

J (interior of E) I X -X' I 

Accordingly, ip = tp^, rather than a definition of the surface, is a supplementary condition 
that must satisfy the solution. 

Although different exact solutions to this problem are known [6] , it can also be posed in 
terms of successive approximations. In the strict sense, this second procedure reproduces 
one of the known exact solutions (the McLaurin ellipsoids). Nevertheless, it is highly 
suitable for the problem we shall study in later sections. 

As is well known, the expansion of (3.12) in spherical harmonics in the region exterior 
to E is: 

«' = -^-E^A(co S 9) (3.13) 



l>2 



where 



Qi= ( R l Pi{cos9) p(R)d 3 R (3.14) 

J (interior of E) 



is the 2'-pole moment. We shall admit that there is a parameter O (O C 1), related 
to the rotation, that measures the deformation of the fluid with respect to the spher- 
ical shape, such that Qi is at least of order Q l . Furthermore, we shall assume that 
it is possible to expand the parametric function of the surface of the fluid, _Re(cos0) 
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in powers of Q. Direct, although tedious, calculation based 
on expressions (3.13) and (3.14) and the equation for E leads to the following results: 
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At the end of the process, it is possible to identify the parameter O, which proves to be: 

RqUJ 2 



o 2 = r2qUj2 



Gm/Ro Gtu/Rq 



(3.16) 



which is the quotient between the rotational energy and the gravitational energy, or, be- 
tween the mean values of the centrifugal force and the gravitational force on the surface. 

4. Relativistic Euler equations. 

Within the framework of relativity, it is also possible to gain information about the 
matching surface from the evolution equations of the fluid. These equations are deduced 
from the energy-momentum tensor conservation condition, which, for a perfect fluid: 



can be written 



T a ? = (p + |r) u a u? + ^g a P (u a u a = -1) 

u a d a p + (p+-|)cU a = 
(p + ^) «%« Q + ^{g aP + u a u p )d pP = 



(4.1) 



(4.2) 



If the metric is stationary and axysimmetric and, also, the velocity of the fluid is in the 
plane spanned by both Killing vectors, it is possible to find coordinates adapted to the 
Killing vectors in which 

u a = </> (5 + ^8%) (4.3) 
and the metric is written in blocks [7]. From the condition u a u a = —1, one deduces that 



- [9tt + 2—gt4> + ~2 9U 
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Under the hypothesis that rotation is rigid (u = constant) and that density and 
pressure do not depend on either t or <p, the second set of Euler equations is reduced to 



d aP = ( P c 2 + p))d a In V {a,b,... = R,S) (4.5) 

Considering this generalisation of the second equation of (3.3) as an equation for the 
pressure, the integr ability condition leads us to p = x(V0> an< ^ the solution is: 

p= / -r^-dX (4.6) 

As in the Newtonian case, p and p are constant over the surfaces ip(R, cosO) = constant, 
and £ (the surface on which pressure vanishes) has as its implicit equation ip(R, cos 6) = 
ipz. If the fluid considered has a constant density, (4.6) leads to: 

P = p(£-l) (4-7) 



5. Parameters of the exterior metric. 

As mentioned in the introduction, the exterior metric that we shall use as the starting 
point is a linear combination of products of factors (2.3) up to R~ 4 . Additionally, the 
interior metric will be calculated by the Post-Minkowskian approximation, considering at 
each order in G a subexpansion in powers of O. To match them to order G, it is of interest 
to know how the terms (2.3) are written as a function of these parameters. Instead of G, 
we shall use the adimensional parameter g = Gm/Roc 2 (Rq is the mean radius of the ball 
of fluid) and we shall consider that the moments have an expansion in O similar to that 
obtained in Newtonian mechanics; that is: 
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We thereby manage to substitute the expansion in powers of R 1 (which lacks sense for 
the interior metric) by an expansion in g and O. This new approach leads to a different 



6 



assignation of orders to the terms of the expansion of the exterior metric. For example, 
according to the first criterion G 2 mB / R A c® and G 2 S 2 /R^c 6 are G 2 / R A terms. However, 
now, whereas the first is of order g 2 fl 2 , the second proves to be of order <? 3 2 . This is a 
direct consequence, in sum, of having chosen O instead of another parameter as an estimate 
of the effect of rotation. Nevertheless, this is the only way to recover the Newtonian results 
at order g. 

Let us now see the meaning of making expansions in these parameters. The parameter 
g is the quotient between the Schwarzschild radius of the body (R s = Gm/c 2 ) and i?o- 
Therefore, g will be small when Rq ^> R s . On the other hand, O has a strange relationship 
with g: 

^I(^) 2 (,3) 

If the body is not very compact (small g), the typical velocity with which a point at the 
surface moves (-Ro^) must be very small as compared with the speed of light, so that a 
small value of Q can be attained and so that the approximation will conserve its validity. 
By contrast, when the body is compact, the surface velocity can be an important fraction 
of the speed of light. In the case of the Earth, the Sun and a typical pulsar (m pu i sar ~ m , 
T ~ 1 s, R ~ 10 km), g and O take the following values: 

Earth Sun pulsar 

g : 1.3 x 10" 9 2.1 x 10" 6 0.14 (5.4) 
O 2 : 1.8 x 10" 3 3.3 x 10" 5 2.9 x 10" 7 

Finally, we can use the expansion in these parameters to determine the matching 
surface in an approximate way from the condition ip(R, cosO) = ips as a function of the 
parameters of the exterior metric (we assume continuity of the metric on E). The final 
result is: 



i? s (cos0) = #0 



1 + (b {2) -\ + \g [5 (0) - l] + 0(g 2 )^j P 2 



+ 0(0 4 ) (5.5) 



It should be noted that since the exterior metric does not contain the octupole moment 
metric (if we assume equatorial symmetry) nor the 16-pole moment, which would go as 
fi 4 , there is little sense in considering expansions beyond O 2 . 



6. The interior metric 

The Post-Minkowskian perturbation method presupposes the existence of a formal 

1/2 

expansion in the gravitational constant G of /i a/3 =( — det(<7 7M )) g af3 — rj al3 [8]. In our 
case, since in the exterior metric there are half-integer powers of g, we shall assume that 
the expansion is in half-integer powers of g (an adimensional parameter). 

n>2 
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Although later, as mentioned above, we subexpand in powers of O each of the terms of 
this expansion, the structure of (6.1) ensures that we can use the whole Post-Newtonian 
scheme, whose details, together with the explicit expressions for the development of the 
different magnitudes appearing in the Einstein equations, can be found in the papers cited 
in reference 7. Accordingly, they will not be referred to here. 

With a view to simplifying the problem as much as possible, we shall assume that the 
perfect fluid is homogeneous (p = constant), in which case pressure will be given by (4.7). 
Since ift ~ 1 + O(g) pressure is order g, which means that the energy-momentum tensor at 
zero-order will be expressed only in terms of p (Newtonian theory). 

Order g. The solution of the Einstein equations at this order, which is regular at the 
origin, is 



with p' = (4tt/3)pRq. The fact that the parameter of the exterior metric m appears in 
(6.2) arises from having chosen g as the constant in which we expand the solution, and has 
no further transcendence. 

Our aim is to match this metric with (2.1) in the sense of Lichnerowicz (see ref. 2): 
the metric and its first derivatives should be continuous on the surface determined by 
(5.5). For the exterior metric (2.1), ^t(i) = ^ext(i) = ^' we tnus cnoose the solution 
of the homogeneous equation mentioned in (6.2) equal to zero. Also, the fact that in the 
exterior metric only terms in Po and P-i appear implies that, by continuity, Hn = if 
n > 2. Since there are four constants, m, B, Hq and H% , and since four conditions are 
imposed by the matching (terms in Pq and Pi in /i9R and their normal derivative to E), 

the matching of the two metrics is possible by identifying the exterior X a and interior x 13 
coordinates. Thus, X a = x a + 0(g 3 ^ 2 ) and the metrics match if: 
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Finally, one has that 
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Order g 3 / 2 . Calculation of the solution at this order passes through the same steps as 
the previous order and in a completely similar way, since only the h^l /2) components are 
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different from zero and, in fact, depend on a single function. Neither is it necessary here 
to modify the interior or exterior coordinates; that is X a = x a + 0(g 2 ). The matching 
conditions demand that 

(6.5) 
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Order g 2 . At this order, the solution obtained 
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solution of the homogeneous equation 
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(where e* = 8\, n l = x l /r y m l = e l j k n J e k ) does not match the part corresponding to the 
exterior metric if one does not make a change of coordinates; this even occurs if O = 
(Schwarzschild). 

The coordinate change 
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induces the change 



h 'a() = h a(3 + g 2 L^q^ - 2d^ a f^] + 0(g b ' 2 ) (6.9) 

in the metric density. In order to preserve the structure and explicit symmetry of the 
metric, one must take: 



f = o 

/• = £// = £ 

l>i i>i 



P} .1 (6.10) 

sin 6 



(Pf are the associated Legendre functions of order 1) where xi and ^i(R) are arbitrary 
functions. The new coordinates X a are still harmonic, for the interior metric, if A/ = 0. 
It is possible to check that the change in the spatial components of the metric density 
that induces the most general solution of this equation, which is regular at the origin 
of coordinates, has the same structure as the solution of the same type and divergence 
free of the homogeneous equation that must be satisfied by the spatial components of 
h^ 2 y Accordingly, we can take this term out of (6.7) since it is included in the change of 
coordinates. 

In order to not introduce into the metric polynomials in cos O of a higher degree than 
necessary, only /| and fl must be taken into account, with the restriction 

3x3 + r 3 = (6.11) 

The ft°2) component, which is not modified by the coordinate change, must be zero in 
order to match the exterior. Also, the matching of implies four conditions and that of 
fo^2) twelve (the term 5^ provides four, n*n J , four, n^e^\ two and e l e J two more). As free 

parameters, one has Hq , H\ , m and B [by the latter two we understand the corrections 
of order g to the Newtonian values (6.3)], which allow us to match with and without 
a coordinate change (see order g), and also the parameters provided by the coordinate 
change. 

The change (6.10) with the restrictions indicated introduces three independent ar- 
bitrary functions xu ^i and X3- If O = (spherical symmetry), (6.9) reduces to f l = 
f(R)N l ; that is, a single arbitrary function that is a particular combination of xi and r±. 
Therefore, all the functions, with the exception of this combination, are at least of order 
O 2 . The constants introduced by the coordinate change are the values of these functions 
and their first and second derivatives on the surface E. At the order of approximation 
required, this is equivalent, for functions of order O 2 , to taking the value of the function 
(and respectively its derivatives) in Ro; for example: 

X3 |e=X 3 ( j R )+O(O 4 ) (6.12) 
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However, for the combination that we have called /, we must conserve the first two terms 
of its expansion in powers of O 



/ | E = f(Ro) - -n 2 f(R )p 2 + o(n 4 ) 
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(6.13) 



and the same for its derivatives. To conclude, we have available 4 + 3 + 3 = 10 arbitrary 
constants (actually, nine because under matching conditions two of them always appear 
in the same combination) to satisfy twelve conditions. As a result, three of them must be 
linear combinations of the rest. This is the case here. 

Having achieved the continuity of the k%s components and their derivatives, and 
therefore having determined the constants associated with the coordinate change, we look 
at the /i°2) component. From the matching we obtain: 



m = p 



1 + 9 



B = p'Rl 



fr 
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(6.14) 



Finally, it should be mentioned that by changing the coordinates not only of the 
interior but also of the exterior metric one can find a global harmonic coordinate system; 
i.e., in these coordinates the metric and its first derivatives are continuous on E. 
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